Multisetting Bell inequalities for N spin-1 systems avoiding KS contradiction 
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Bell's theorem for systems more complicated than two qubits faces a hidden, as yet undiscussed, 
problem. One of the methods to derive Bell's inequalities is to assume existence of joint probability 
distribution for measurement results for all settings in the given experiment. However for spin-1 
systems, one faces the problem that eigenvalues of observables do not allow a consistent algebra if 
one allows all possible settings on each side (Bell 1966 contradiction), or some specific sets (leading to 
a Kochen-Specker 1967 contradiction) . We show here that by choosing special set of settings which 
never lead to inconsistent algebra of eigenvalues, one can still derive multisetting Bell inequalities, 
and that they are robustly violated. Violation factors increase with the number of subsystems. The 
inequalities involve only spin observables, we do not allow all possible qutrit observables, still the 
violations are strong. 



Introduction: Intrinsic randomness of Quantum Me- 
chanics has always been difficult to accept. It was chal- 
lenged by Einstein, Podolsky, and Rosen [1] . Many years 
later Bell [2] formulated an inequality invalidating this 
challenge. Violations of Bell-type inequalities in an ex- 
periment or theory exclude local and realistic comple- 
tions of quantum mechanics. There have been many ex- 
perimental attempts to violate a Bell inequality, for a 
recent review see [3], often related to some specific chal- 
lenges. For example, in Refs. jlHS] the challenge was 
to close various loopholes allowing a Local Realistic de- 
scription. In the meantime, it turned out that Bell's 
inequalities are useful tools to classify non-classicality of 
quantum correlations, which could be applied in vari- 
ous quantum information protocols breaching the lim- 
its of classical communication and classical computation. 
Thus, Bell inequalities have now a practical importance. 

We have many generalizations of the inequalities. For 
more than two qubits Mermin introduced a series of 
inequalities for experiments with choice between two 
dichotomic observables, [7]. Alternative versions were 
given in Refs [El |9] , and finally, the full set in [10l HI] , 
and [12] . More recently, a further generalization, to more 
than two settings per observer, was introduced e.g. in 
Refs. 13J and [H]. 

In the case of higher-dimensional systems early results 
can be found in Refs [15]. Later, Kaszlikowski et al. [T5] 
have provided an evidence that two entangled qutrits vio- 
late Bell inequalities stronger than two entangled qubits. 
This was confirmed in [18 . A Clauser-Horne-type [17] 
inequality was proposed for two qutrits in Ref. [19]. Fi- 
nally, Collins et al. have derived series of two-setting 
two-observer Bell inequalities involving measurements of 
any finite dimensionality |20j . and a different form has 
been derived by Zohren and Gill, [21 . A numerical study 
of more settings scenarios two-observer scenarios can be 
found in Ref. [22]. 

A three-dimensional Hilbcrt space is sufficient to con- 
duct a proof of the Kochen-Specker (KS) theorem [23] . 
which falsifies non-contextual hidden variables. Recently, 



a non-contextual realistic description of a single three- 
dimensional system was falsified in an experiment |25j . 
This may make constructing Bell inequalities for qutrits 
look pointless. Nevertheless, a closer examination of the 
problem still provides us with a motivation. The KS 
theorem shows us that there cannot be hidden variables 
describing results of all measurements, as some subsets 
of measurement cannot be modeled so. Such sets, of- 
ten called KS sets can be now found algorithmically [24] . 
However, there are certain other subsets, for which these 
models are not invalidated by KS-type reasonings. 

Here we introduce multisetting Bell-type inequalities 
for collections of spin-1 systems, which we shall call cone 
Bell inequalities. We shall consider only observables re- 
lated with spin measurements, this is a completely dif- 
ferent approach than in [TBI H5H2D] . where observables 
with no classical analogues were considered (the specific 
observables involved multiport beamsplitters). This is 
the reason why we call our subsystems spin-1 systems. 
The terminology stresses that we do not allow all possible 
qutrit observables. The word "cone" used to describe the 
inequalities reflects the fact that the Stern-Gerlach-type 
measurements entering the inequalities have the prop- 
erty that the vectors defining the measured component 
of the spin can come only from a specific cone (described 
in more detail further on). The cone sets are such that 
every projective measurement appears in a single (local) 
context, thus they cannot form KS sets leading to a con- 
tradiction, for whatever the number of settings. One can 
fully separate the two no-go statements against hidden 
variables. And perhaps even more important is the fact 
that the shown Bell-type violations of local realism are 
robust. 

Conceptually, these inequalities are related to earlier 
works on qubits 26-29 . An exponentially strong viola- 
tion of the inequalities will be demonstrated. 

We will utilize only measurements, for which no KS 
argument can be constructed. Note that therefore our 
results apply also to the case of locally contextual hidden 
variable theories. 
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Observables: Our construction will rely on squared 
spin-1 components (S ■ ft) 2 , where n is a unit vector. 
These observables are particularly interesting for their 
association to the sphere of Stern-Gerlach "magnet" ori- 
entations. It should be stressed that this sphere is differ- 
ent from the Bloch sphere parameterizing qubit observ- 
ables. In the case of squared components of spin-1 we 
have identification n f-> —ft. 

Since the square of the total angular momentum of a 
qutrit is 2, we have (we put h~ = 1) 

S 2 = (S ■ x) 2 + (S ■ y) 2 + (S ■ z) 2 = 2. (1) 

The three spin components mutually commute and thus 
obey the 1-0-1 rule, which gives rise to all KS argu- 
ments for qutrits. That is, a simultaneous measurement 
of all three squared orthogonal spin-1 components will 
always give two l's and one in some order. 

Note that one can introduce state |n), for which (S ■ 
fi) 2 \n) =0. A kind of Malus law for spin-1 systems reads 

|(n|n')| 2 = (ri-n') 2 . (2) 

The whole set of squared spin-1 components is still 
too large for our purposes. There cannot exist hidden 
variables assigning outcomes to all measurements on the 
sphere. A simple proof of this statement was described in 
Ref. [30]. Bell's theorem would formally not even make a 
sense. Hence we focus only on observables parameterized 
by vectors n(<p) — ^= ( V2 cos 0, v2sin0, 1). In such a 

case (S ■ n((f>)) 2 is represented by a matrix 

(2 er 1 ^ e~ 2i ^ \ 
e«* 2 -e-** . (3) 
e 2# _ e i4> 2 / 

Note that any three observables S ■ n(4>j ) with cj>j — <fio + 
jy, where j — 0, 1, 2 constitute a triad of measurements 
at mutually orthogonal directions and hence their squares 
obey the 1-0-1 rule. Moreover, two different such trios, 
with <po — (f>' ^ j™, do not share any eigenvectors and 
hence are completely independent from each other in a 
local realist's eyes. 

Construction: We shall now demonstrate how the con- 
flict between the local realistic and the quantum mechan- 
ical descriptions of entangled qutrits may arise. The 
derivation is most concise if one uses instead of opera- 
tors S 2 (4>) traceless observables, 0(4>) — S 2 (4>) - 2/3. 
Obviously, their eigenvalues are 1/3 and —2/3. 

In every run of an experiment let each of observers 
1,2,3,. ..,7V receive one spin-1 system from a common 
source. Each of them measures a single squared spin com- 
ponent in direction n(<pi),n(<p2),n(4>3), respectively 
(getting | or — | , and thus obtains one of the eigenvalues 
of his/her observable 0(<fr). Subsequently they average 
the product of their results over many runs. This aver- 
age will be called the (quantum mechanical) correlation 



function, Eqm{<I>i, <f>2, 03, ■•■) and is an element of some 
real linear space of functions. In this space we introduce 
a standard scalar product, 

{A,B)= M<h,<h,~)B(<h,<h,~) (4) 

4>i € Si 

<t>2 G S 2 

in the discrete case (Si, S2,... are countable), or 

{A,B) = J 0ieSi #i4 eS2 #2-. 

xA(&, &,...) (5) 

in the continuous case (Si, S2,... are not countable). 
Hereafter, we will rather use the sum notation, but all 
considerations are valid for both cases. Here we focus 
only on the symmetric case of Si = S2 = ... = S. To 
take a benefit of the 1-0-1 rule, we assume that whenever 
some (f>o belongs to S, so do (f> + ^ and 0o + 

The usual assumptions of local hidden variable theories 
lead to following structure for correlation functions 

N 

£lhW>i,</>2,...)= dXp^HhifaX), (6) 

J 1=1 

where A represents hidden variables, p(X) is their distri- 
bution, and each Ii(4>i) takes values equal to the local 
eigenvalues, that is | or — |. The "unusual" additional 
assumption is that Ii(</>i) + 1% + ™) + h [cf>i + ^p) = 
0, which expresses the 1-0-1 rule. Such models form a 
convex set e^n (the additional assumption plays no role 
in the case of this property). For its extreme points 
all measurement outcomes are deterministic, and the 
whole correlation function £?iij(<^i, <p2, •■•) factorizes into 
Yli—i Ii(<j>i)- If the following inequality holds for all 
Elr € €lb- 

e qm — ( E QM,Eq M ) > (Eq M ,E lr ), (7) 

then Eqm 7^ Elr, that is a local realistic description 
is not possible. Notice that since the right-hand side 
of inequality ^ is linear in Elr, its maximum will be 
attained for one of the extreme points. 

A simpler form of the inequality: One can simplify 
the expression for max£ tJi (Eqm, Elr)- Let us consider 
the example of 2 spin-1 systems. Since Ii(4>i) takes values 
I and -§, we can write = \ — Xi{4>i)i where x»(<&) 

is a characteristic function giving values 1, if 4>i belongs 
to the subset for which Ii(4>i) = — I, and otherwise 0. 
Thus e.g., 

ma,XE LR (Eqm, E L r) 

= E^esE^es^QMO/Ji^X! -Xi(<M)(! ~ X2(h)) 

(8) 
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However, 



Ea M {<t>i,4>2) 



y 

2 es 



^ EQM{<t>i><h) = 0- 



Simply from the 1-0-1 rule any trio of the 
form EQ M {(t>i,4>2), E QM (0i + § tt, fo) and 
Eqm {<t>i + f 71 ") 02 ) sums up to 0. Thus the only 
term that does not vanish is the one with the product of 
the characteristic functions. Therefore, we get 



ma,x(E QM ,E LR ) - , , 

Elr 



E 



E QM {<t>i,fo), (9) 

where o~i is the set of values of fa, for which Ii{<fii) = — |. 
All this generalizes to the multiparticle case in a trivial 
way. 

Falsification of local realism: We will demonstrate in- 
stances of violation of inequalities of the type ^ . Let us 
focus on biased GHZ states, 



IV 



N 



+ ±10)®^ + ) , (10) 



where | ± 1) and |0) are eigenstates of S z with respective 
eigenvalues ±1 and 0. The quantum mechanical correla- 



tion function for state (10) and observables 0(0) reads 



3 2 + ' 



£qm(01,02,-) 



(ii) 



The values of {Eqm, Eqm) for the continuous case, 
£ = [0, 2tt) read 



{Eqm , E, 



QM) 



64(2tt 



3 (4+JV) 



(12) 



whereas for a uniform distribution of 3n settings, S 



2j7T 

3n 



J 



0, 1, 2, ..3n — l} we have 



{Eqm, Eqm) 



64rt 



N 



3 (4+JV) 



(13) 



(the formula does not hold for the irrelevant case of n = 
!)• 

Estimation of upper bound for local realistic models, 
for specific sets of settings: Here we are going to present 
the method to bound from above {Eq M ,E lr ). First, let 
us introduce 



- 2lT k 2n ' 
3n 3 ^ 



(14) 



as settings that we shall work with. Here k = 1, ...,n, 
I = 0,1,2, whereas x — 1,2, ,..,N denotes the observer. 
We want to estimate the maximal value of 

{Eqm, E lr ) 



E{^ j \...,^ j nx,M lj nXaM 2h )---X 



Xo-xifix) i s the characteristic function, equal to 1 when 
G Gi and otherwise. Note that for any k there is 
always only one j for which Xa !C {4'x) = 1- This is a 
manifestation of the 1-0-1 rule. 

We can decompose the correlation function according 
to well-known expressions for trigonometric functions. 

Eqm^i, ■-, 0Jv) 

= 3 2% N (COS 0i COS 02 COS 03 . . . COS ■ ■ ■ 

-K-l^cos 20i cos 20 2 ... cos 2(f> N + ...)), (16) 



where 



COS Ipl COS 1p2 COS COS 1p N 



(17) 



stands for the sequence forming the decomposition of 
cos^Vi- As we shall see, the actual details of such 
a decomposition will be irrelevant. What is important is 
that we have a certain sequence of terms in form of prod- 
uct of trigonometric functions, each one containing either 
cosipk or sinipk- Each term is multiplied by a coefficient 
±1. 

Notice that the values of cos tj^P , sin 0^p cos 20^p , and 
sin 24>x form four orthogonal 3rt-dimensional vectors c%, 
s*i, C2, and s*2, respectively, which all have the same equal 
norm, denoted here as M . Hence we can put problem 



(15) into the following form 



m& *E LR { E QM,E LR ) 

max (Tli .. (Tiv (-l) Ar [ci <g ci <g .. 
+ {-l) N {c 2 ® c 2 (g ... ® c 2 + ...)] • xi 



32+N 



C\ + ... 



where Xx is a 3n dimensional vector with entries 0, 1 cor- 
responding to the characteristic function XaS 1 ^), an d 
• represents the standard real "dot" scalar product in 
a 3nN dimensional real space. The compound symbol 
ci®Ci®... ®ci + ... represents a sequence of tensor products 
in obvious relation with the former sequence of trigono- 
metric functions (17). That is every cos is replaced by 
ci , every sin by si , and the coefficients of the expansion 
are kept unchanged. 

Now we use an important technical observations. Fist 
we notice that, 

(q ® q (g ... (g) ci + ...) • xi <g ■■■ <g Xn 
= (Q • Xi){ci ■ x.2)--{ci ■ xn) + 

(18) 



where now • stands for a scalar product in a real 3n di- 
mensional space. Since the four vectors, ci, si, C2, and 
S2, are orthogonal, we can always put 

Xx-ci = ||xjj||Mcos^cos(a x + Z§) , (19) 

Xx • = \\j&\\M an?* caa (a, + if) , (20) 

where / = 1,2 and a x and are some angles, and 
(15) \\xx\\ denotes the norm of a projection of x. x onto the 
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4-dimensional space spanned by c\, s±, <?2, and s*2- This 
is an obvious property of any decomposition of a projec- 
tion onto a 4 dimensional subspace. If whatever vector x 
is projected into a subspace spanned by normalized ba- 
sis vectors oi, b\, et2, &2, one can always find such three 
angles j3±, /?2 and a so that one can put 



X' 



2 



|x" || cos(a ■ 
1=1 



\l)(cos foai +smj3ih). (21) 



After putting the above formulas to ( 18 ) we get the 
following 

(E QM ,E LR )<^r(M\\£\\) N 
x max [nili cosa i C0S (XwLi Z 3 ]) 



(-l) w n 4 =iCos(a l + f)cos(E; V =1 /3: 



JV 



(22) 



The functions of Yli=i /^i result from a reverse ap- 
plication of the expansion formula for trigonometric 
functions of many angles: cos j3[ cos /3 2 cos j3 3 ... cos j3 l N + 
... = cosQ^iLi Pi)- Since cosine is never greater than 
1, the expression in the square bracket can be ma- 

which in 



jorized by 



rJY 



turn can be seen as a scalar product of two vectors, 
(± cos ct\, ± sinai) and ^11^=2 cos > 11^2 sm J > 
both with norm not exceeding 1. By the Cauchy in- 
equality, the value of this expression cannot exceed 1. 
therefore we get 



(EqM: E LR ) < 



32+jv 



(M\\m\) 



(23) 



Thus once we know for the maximal possible value of 
||Xa;||Af, we can estimate the local realistic bound from 
above (note that the above bound does not have to be 
tight). 

For n = 3 we found it by a MATHEMATICA case-by- 
case study to be 



max (Ti ||x[ l ||M=^(l + 2cosf) 
« 2.86822. 



(l + 2sinf|) 



(24) 



Hence this particular estimate can be used for an arbi- 
trary even N. The strength of violation can be then 



estimated from Eq. ( 13 ) as 



< 



{Eqm ,Eq m ) 
inaxEr r (Eqm ,E La ) 



9 I V / ( 1 + 2cOS Tf) 2 + ( 1 + 2sin Tg) 2 



N 



(25) 



which predicts violation of the inequality for N > 3. We 
have found the values of max CTi | \xj \\M for n = 4, 5, 6, 7 



FIG. 1: The violation ratio ^ in function of the number of 
local measurement trios n for 2 spin-1 systems. In all figures 
the solid line shows the continuous limit. In the shaded region 
the results are calculated explicitly, outside this region we 
show an estimate based on the conjecture (see text for details) . 



FIG. 2: The violation ratio -k in function of the number of 
local measurement trios n for 3 spin-1 systems. 



to be 3.7678, 4.678, 5.5932, and 6.5112, respectively. 
Note that these values are smaller than n, which when 



one compares ( 23 ) and ( 13 ) guaranties falsification of Lo- 
cal Realism for sufficiently large N. Interestingly, for 
N = 2 the bound given by Incq. (23 1 is tight for all n. 
Numerical analysis: 

For N — 2 we have numerically tested all local realis- 
tic models for n < 7 and verified that the highest values 
of (E QM ,E LR ) for n = 2,3,4,5,6 are 0.1975, 0.6083, 
1.1852, 1.6632, and 2.380, which corresponds to viola- 
tion ratios h = {EQA ff QM l , equal to 0.8889, .9724, 

1.0018, 1.016, and 1.023, respectively. V tells us that 
to what extend we can maintain the violation, while we 
mix the state IVO^I with the white noise, p no ise (propor- 
tional to the unit operator) , so that we get a mixed state 
P\ip)(ip\ + (1 — P))p n oise- In such a case P = V is the 
threshold value at which the inequalities are not violated 
anymore. Please note that the correlation functions for 
the noise state vanish. 

Similarly we processed all local realistic models for 
N = 3 and n = 2, 3, 4, 5. We have found that i equals 
0.889, 1.401, 1,398 and 1.397 respectively. The results are 
presented in Fig 2. 

Conjecture: All these numerically analyzed instances 
reveal a specific pattern of the form of the optimal sets a 
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FIG. 



3: The violation ratio h in function of the number of 



local measurement trios n for 4 spin-1 systems. 




FIG. 4: The violation ratio h in function of the number of 
local measurement trios n for 5 spins- 1 systems. 



because three of four of these points are associated with 
<Tj's of the second type. Calculations show that the plot 
of the the ratios of (Eqm , Eqm) to maximal local real- 
istic values under the restrictions that oVs are of these 
type has some kinks for mod(n,4) = 1. In Fig. 4 the 
third point corresponds to first such kink, but this value 
is overtaken with a value for cr^'s of the first type. 

These observations we can be extended to the contin- 
uous set of observables. One can take a,; 
(i 7^ N) and &N - 



2tt 
3 



27T 27r 1 

3 ' 3 J 



(Eqm, E lr ) = f_ 3 ^ <%... jjL d^Ar-i f_ 3 2?l„ ' 

3 3 3 ' " 

X g#F2 (cos J^j 4>j - COS 2 J^j 4>j) 
cos(n-(^-l)7r) cos2(a-(^-l)7r) 



37 



3T 



while for even N we have 

2tt 2tt 2tt i a 

{E QM , E lr ) = fj>^ JjW d(j> N _i J_ 3 2^ +a ' 



3W + 2 



cos(a-^Tr) cos2(a-^Tr) 



37 



-+ 2 9iV-3 



3/V 



(26) 



(27) 



(that is those that lead to the highest value of ([9])). There 
are always many local realistic models maximizing the 
scalar product with the quantum mechanical function, 
but by the symmetry of Eqm (which is solely a function 
of X)i=i 4>i) one °f them has <7i = a% = ... = ctn-i = & 
and ct/v = R a {°~)- Ra is an operator, which shifts its 
argument by a modulo 2tt, R a (b) = a + b. Two kinds 
of ers are particularly relevant. One, denoted as a' , has 



-4-2/oiV 



(2' 



[0. 



or a 



all points bunched together, e.g. a 
{0, i^}. The other, called a", has two equally 

sized parts with possibly equally sized gaps, e.g., a = 



[of)uhf) 



{o 



2tt 
9 



4rr IOtt 
9 



12tt 14tt 



In all explicitly studied cases the optimal sets <7j for 
the local realistic value were of types either a or a" . This 
observation can be taken as a conjecture for more compli- 
cated situations. For N — 2, 3 we have been able to give 
estimates of maximal violation under this assumption to 
up to 30 trios of observables. This can be seen in Figs. 1 
and 2 and as points outside the shaded region 

For higher N we mainly used oVs of either a' or a" 
type. The results are presented in Figs. 3 and 4. Wc 
have been able to estimate violations for N < 6. This 
hypothesis was additionally supported by studying all 
local realistic models for N = 4, n = 2,3,4 and N — 
5,n = 2,3. 

Notice that the points in the plot for N — 5 are dis- 
tributed in a different way than in other plots. This is 



The maximum of these functions is 2 3 N 3 
1). This, when compared with (12) corresponds to a 

8 



violation ratio -k 



However, this is 



9(2« + l) V3V37 
only a conjecture. 

Closing remarks: We have proposed a new family of 
Bell-type inequalities for arbitrarily many spin-1 systems. 
Observers might have a divisible by three, but other- 
wise an arbitrary number of observables per side. These 
inequalities totally avoid the Kochen-Specker contradic- 
tion, but still use the 1 — 0—1 rule, which is the source of 
it. Despite specific constraints that we introduced to our 
analysis, violations of these inequalities are quite high, 
and grow exponentially with the number of spin systems. 

On the trail of the derivation of the inequality we have 
tried to introduce some free parameters. For example, 
we have given up the tracelessness of the observables. 
However this only decreased the strength of violation. 
Similarly, changing the the degree of the biasedness of the 
state (i. e., the relative weights of its three components) 
had only a negative effect on the strength of violation. 
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